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Abstract. During division it is of primary importance for a cell to correctly 
determine the site of cleavage. The bacterium Escherichia coli divides in the center, 
producing two daughter cells of equal size. Selection of the center as the correct division 
site is in part achieved by the Min-proteins. They oscillate between the two cell poles 
and thereby prevent division at these locations. Here, a phenomenological description 
for these oscillations is presented, where lateral interactions between proteins on the 
cell membrane play a key role. Solutions to the dynamic equations are compared to 
experimental findings. In particular, the temporal period of the oscillations is measured 
as a function of the cell length and found to be compatible with the theoretical 
prediction. 
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1. Introduction 

During division a cell has to specify in particular the location of cleavage. In the rod- 
shaped bacterium Escherichia coli, the division plane is determined by the location of the 
Z-ring p. This structure is built from FtsZ-filaments and forms on the inner bacterial 
membrane. The position of the Z-ring in turn is first of all determined by the distribution 
of the nuclear material inside the cell. A mechanism termed "nucleoid occlusion" 
restricts the formation of the ring to regions void of DNA [2] Ej- After duplication 
and segregation of the chromosome, three locations of possible ring formation remain: 
at the cell center and close to the two cell poles. Selection of the center as the correct 
division site is achieved by the Min system 0] 13] . Deletion of any of the Min proteins 
results in division septa forming close to one of the two cell poles in about 50% of all 
divisions. In these cases, DNA- free mini-cells are formed jHj- 

The Min system consists of three proteins, MinC, MinD, and MinE. Out of these, 
MinC induces the depolymerization of FtsZ-filaments and inhibits the formation of the 
Z-ring [7j. The distribution of MinC on the membrane changes periodically with time 
such that in one half of the cycle, MinC accumulates at one pole while it accumulates 
at the opposite pole in the second half of the cycle [HI E]- Formation of the Z-ring 
is thereby suppressed at the cell poles. The temporal period of the oscillation ranges 
between 40s and 120s in wild-type cells. In bacteria of a length that exceeds a certain 
threshold, a striped oscillatory pattern appears, where the number of stripes increases 
with increasing cell length. This observation is indicative of an intrinsic spatial wave- 
length of the oscillations. 

The oscillations of MinC require the presence of both MinD and MinE, which 
themselves also oscillate ^01 ITT] . In fact, MinC binds to MinD and follows its 
dynamics [TO]. In contrast, MinE is mostly localized in a ring structure which oscillates 
around the center of the bacterium. Remarkably, MinC is not necessary to generate 
oscillations, as MinD and MinE oscillate also in the absence of MinC. The behavior of 
MinD and MinE has been characterized by intensive biochemical and genetic studies. 
In vitro experiments have shown that the ATPase MinD has a high affinity for the inner 
bacterial membrane if ATP is present ^2]- For concentrations of MinD exceeding a 
critical value, filamentous MinD aggregates are formed on the membrane [T2~] HHj. The 
formation of MinD aggregates is likely to be a two step-process, where MinD first binds 
to the membrane and then self-assembles [12] . Indeed, in presence of ATP7S, a non- 
hydrolysable analog of ATP, MinD associated with the membrane but failed to form 
filaments. As for MinE, it associates with the membrane only in the presence of MinD. 
There, it stimulates hydrolysis of the ATP bound to MinD, which eventually drives the 
proteins off the membrane |12j . 

These in vitro results are compatible with the behavior of MinD and MinE in 
vivo. In MinD depleted cells, it was observed that MinE is dispersed in the cytosol, 
while MinD is homogenously distributed on the cytosplasmic membrane if MinE is 
absent jlOj . Furthermore, helical MinD aggregates have been observed to form on 
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the inner membrane [Hj. The significance of the helical structures for the oscillation 
mechanism is still not understood. Finally, the oscillations do not depend on the 
synthesis and degradation of the Min-proteins [HI] . 

Theoretical investigations of the Min-system suggest that the periodic transloca- 
tions of the Min proteins can be attributed to a collective effect of many interacting 
molecules resulting from a dynamic instability [13 dl 113 CHI • Central to all proposed 
mechanisms is the attachment of MinD to the cytoplasmic membrane, recruitment of 
MinE to the membrane by MinD, and dissociation of MinD from the membrane induced 
by MinE. The mechanism proposed by Meinhardt and deBoer J3| belongs to the class 
of classical reaction-diffusion systems with short-range activation and long-range inhibi- 
tion. The synthesis and degradation of the Min proteins play an essential role. Howard 
and colleagues ^Hj assume that MinD and MinE form complexes in the cytoplasm, 
which then bind to the membrane. Membrane-binding is hampered by MinE present on 
the membrane. Furthermore, the protein number is conserved. The same holds for the 
mechanism presented in Ref. [Tj\. There, however, first MinD binds to the membrane 
and then recruits MinE. More importantly, aggregation of membrane-bound MinD is 
essential. In contrast to react ion- diffusion systems, the instability is here driven by the 
aggregation current of MinD. Rather similar to this mechanism is the one proposed by 
Huang et al. ^H] more recently. In contrast to Ref. jTTj, aggregation is there assumed 
to be a consequence of MinD binding cooperatively to the membrane. This seemingly 
small difference in the formation of membrane-bound MinD aggregates has remarkable 
consequences. Firstly, it is essential to describe the Min dynamics in a three-dimensional 
geometry. Secondly, a finite ADP to ATP exchange rate for cytosolic MinD is a key in- 
gredient. As transport is purely diffusive, the instability leading to the oscillations is in 
this case essentially of the same kind as in the other reaction diffusion systems [TBI I16j . 

In this work we re-investigate the mechanism proposed in Ref. J7j. There, the 
aggregation of membrane-bound MinD was formulated in terms of a kinetic hopping 
model. Here, we will use a phenomenological description, which allows for a quantitative 
comparison with experimental results. The paper is organized as follows. First, we will 
describe the equations governing the dynamics of the protein distributions in the cytosol 
and on the membrane. We then analyze the system in the limiting case of homogenous 
cytosolic protein distributions and discuss the oscillatory solutions. The dependence of 
the temporal oscillation period on the system length is compared to experimental data. 
Afterwards we discuss possible mechanisms underlying the formation of the MinE-ring. 
Finally, we discuss our results in relation to the other proposed mechanisms as well as 
implications for possible future experiments. 

2. Dynamic equations 

As mentioned above, the periodic changes in the distributions of the Min proteins require 
the presence of MinD and MinE, but not of MinC. Therefore, we will focus in the 
following on the dynamics of MinD and MinE. Motivated by the observations reported 
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Figure 1. Schematic representation of the dynamics of MinD and MinE. a) 
Exchange of MinD and MinE between the cytosol and the membrane. 1) MinD 
undergoes a conformational change upon binding ATP, 2) ATP-bound MinD binds 
to the membrane, 3) MinE binds to membrane-bound MinD, and 4) MinE-induced 
ATP hydrolysis leads to detachment of MinDE-complexes from the membrane, b) 
Interaction of membrane-bound proteins leads to the formation of MinD aggregates. 



above, the dynamics of the Min proteins is assumed to be driven by four properties 
of the Min proteins i) a high affinity of ATP-bound MinD for the membrane, 

ii) a high affinity of MinE for membrane-bound MinD, iii) a MinE-induced increase of 
the ATP hydrolysis-rate by MinD, which leads to the detachment of MinDE-complexes 
from the membrane, and iv) interactions between membrane-bound proteins. The last 
property accounts for the formation of MinD aggregates on the membrane, which is 
likely to result from self-assembly of membrane-bound MinD ^2] • In addition, proteins 
are transported by diffusion. A schematic representation of the Min dynamics is given 
in Fig. □ 

Formally, the dynamics is given in terms of the concentrations of cytosolic MinD 
and MinE, cd and ce\\, as well as the concentrations of membrane-bound MinD and 
MinDE-complexes, Cd and Cd e - In the direction perpendicular to the long axis of the 
bacterium, diffusion homogenizes the cytosolic distributions on time scales that are 
short as compared to the temporal oscillation period. Assuming in addition that MinD 
aggregates into a linear structure on the membrane, the dynamical equations for the 
protein densities in the cell can thus be reduced such that they depend only on the 
position x along the long axis of the bacterium, see | Appendix A| Explicitly, 

9 t C D = - CUo(c max — Cd — C de )c D + UJ de Cde + D D d 2 x C D (1) 

d t c E = u de Cde ~ u E c d c E + D E d 2 x c E (2) 
d t c d = uJ D (c mSLX - c d - c de )c D - uj E c d c E - d x j d (3) 

d t C de = - U de Cde + U E C d C E - d x j de (4) 

The properties (i)-(iii) lead to an exchange of MinD and MinE between the cytosol and 
the membrane. The corresponding reactions are described as first and second order 

|| MinE forms dimers an d ce is actually the distribution of MinE dimers. In the following, the 
term "MinE molecules" refers to these dimers. 
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processes. The density of available binding sites for MinD on the membrane is given by 
Cmax — Cd~ c de, where c max is the maximal possible value for the protein density on the 
membrane, and ^(cmax — Q — Q e ) is the binding rate of MinD to the membrane. The 
binding rate of MinE to membrane-bound MinD is o^eQ, while Ud e is the detachment 
rate of MinDE complexes from the membrane. We assume the complexes to consist of 
one MinD and one MinE molecule. Dp and De are the respective diffusion constants 
for cytosolic MinD and MinE, and the interactions of membrane-bound proteins are 
captured by the currents jd and jd e - Note, that in these equations the rebinding of ATP 
to MinD after detachment from the membrane is assumed to occur on a sufficiently 
short time-scale such that it does not need to be incorporated explicitly. The effect of 
a finite ATP exchange rate will be discussed below. 

The current of membrane-bound MinD has a diffusive part and a part due to 
the interaction between MinD proteins. In order to capture generic effects of the 
interaction, the current of membrane-bound MinD is taken to be of the Cahn-Hilliard 
form. Explicitly, 

jd = -DddxCd+Cdic^-Cd-Cde^hdxCd+hdlcd + hdxCde + hdlcde}- (5) 

In this expression, Dd is the diffusion constant of the MinD proteins on the membrane 
and the coefficients k± and &2 are phenomenological parameters that describe the 
interaction between MinD molecules. Possible modifications of this interaction due to 
the presence of MinE are taken into account by the parameters k\ and ki that describe 
the interaction between MinD and MinDE-complexes. Note, that for an attractive 
interaction k\ > 0, while k\ < in the opposite case. Stability on small length scales 
requires &2 > 0. The current of MinDE complexes has the same form, but for simplicity 
will be omitted in the following. 

Finally, the boundary conditions have to be specified. We impose zero flux at the 
boundaries, such that the total protein numbers 

rL/2 

/ dx (c D + c d + c de ) = LV (6) 

J -L/2 
rL/2 

/ dx (c E + C de ) = L8 (7) 

J-L/2 

are conserved. Here, L denotes the length of the system and LT> and LS are the total 
numbers of MinD and MinE molecules in the system, respectively. 

3. Homogenous cytosolic distributions 

We now analyse the dynamic equations (|Tj) -P| in the limiting case of homogenous 
cytosolic MinD and MinE distributions, i.e., cr>(x, t) = cdf) and Ce(x, t) = c^(t). This 
corresponds to the case where the times needed for MinD and MinE to diffuse along the 
whole length of the bacterium, L 2 /Dd and L 2 /De, respectively, are short as compared 
to all other relevant time-scales involved. In this case, the dynamics of the cytosolic 
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distributions is described by ordinary differential equations 

^c D = - w D (c max - V + c D )c D + 0J de (£ - c E ) (8) 

^-c E = - uj e (D - £ - c D + c E )c E + u de (£ - c E ) . (9) 

Here, the distributions of membrane-bound MinD and MinDE have been eliminated 
using equations © and (JJJ). 

Under the conditions < Cd < T> and < c E < £ , the above equations have one 
and only one fixed point. This point is always stable and, asymptotically, the cytosolic 
distributions will approach the corresponding stationary values C E and C E , respectively. 
In this limit, the dynamics of the Min proteins is described by two partial differential 
equations for the distributions of the proteins bound to the membrane 

d t c d = uJ D C D {c max -c d - c de ) - u E C E c d - d x j d (10) 
d t c de = - u de c de + u E C E c d . (11) 

Note, that the reaction terms in these equations are linear and describe relaxation 
to a stationary value; only the current contains non-linearities and can generate an 
instability. This feature distinguishes this system from classical reaction-diffusion 
systems, where transport is due to diffusion and where instabilities are created by the 
reaction terms. 

The homogenous state c d (x) = V — £ — Co + C E and c de (x) = £ — C E is a stationary 
state of the dynamic equations (fTU|) and (JTTJ). It is stable, unless ki exceeds a critical 
value fci )C . The results of a linear stability analysis for a supercritical value of k\ are 
shown in Fig. |2K- The stability region of the homogenous state as a function of the 
total MinD and MinE concentrations, T> and £, is shown in Fig. |2b- At the instability 
an inhomogenous stationary state appears if the detachment rate of MinDE complexes 
from the membrane is above a certain critical value, uo de > uj de>c . In the opposite case, 
an oscillatory state appears. Oscillatory instabilities only occur if the protein density 
on the membrane cannot exceed a maximal value c max . For an oscillatory instability the 
unstable mode is of the form 

c d oc cos(Q c t) cos(g c x) (12) 
c de oc cos(f2 c t + (p) cos(g c x) (13) 

This standing wave reflects the qualitative features of the observed Min-oscillations. 
The wave number q c = rnr/L, where n is a natural number, and the frequency Q c of the 
critical mode depend on the system parameters. For instance, we find 

^4 _ {upC E + U de + U) E C E ) 
Cd(Cmax — C d — C de )k,2 

and if ki = /c 2 = 

Q 2 C = u D u E C D C E - u 2 de (15) 

The oscillatory patterns can be obtained from numerical integration of the dynamic 
equations (|1U|) and (|llj). A typical example is shown in Fig. |2t, b. For some time the 
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a) b) 




Figure 2. Linear stability of the homogenous state, a) Real (Re, solid line) and 
imaginary part (Im, dashed line) of the eigenvalues of the linear operator describing 
the dynamics of small perturbations around the homogenous state as a function of 
the wave number q. Modes with wave numbers between l/im -1 and 2.2/zm -1 are 
oscillatory and unstable, b) Stability of the homogenous state as a function of the 
average total MinD and MinE densities V and £. The solid line indicates a line of 
oscillatory instabilities while the dashed lines indicate stationary instabilities. The 
values of the parameters are loo = 4 • lO^/ims" 1 , loe = 3 • 10 _4 /ims _1 , u>de — 0.04s" 1 , 
D d = 0.06/rai 2 s -1 , c max = 1000/im- 1 , ki = 1.5 • lO^Vm 4 * -1 , k 2 = 1.8 • 10" 7 /im 6 s" 1 , 
k x = -1.2 • 10- 6 /im 4 s~\ h = 1.2 • 10- 10 /xm 6 s- 1 . In (a) V = 900^-* and 
£ = 350^m- 1 . 

total MinD-distribution q + Cd e is localized in one half and then switches to the other. 
In this process, the transition time is very short as compared to the dwell time in one 
half. The MinE distribution shows a similar behavior, but the transition between the 
two halves is less rapid. The time-averaged distribution of both, MinD and MinE shows 
a minimum in the center and increases towards the system boundaries, see Fig. Efc. The 
parameters have been chosen such that the temporal period is about 80s, which is similar 
to the values observed in experiments with fluorescently labeled MinD, see Fig. HJ The 
figure also displays the time-averaged MinD-distribution with a minimum in the center. 
In the case displayed on Fig. Ef, the minimum at the center is more pronounced than 
for the theoretical calculation: while experimentally the minimum is at about 50% of 
the maximum, it is at about 70% in the numerics. This might indicate the need for 
further non-linearities in the theory. However, for other cells examined, the minimum 
is much shallower or even absent (data not shown). This might reflect deviations in the 
total protein density in individual bacteria from the average total protein density in a 
bacterial colony. Note also, that in the numerics, the value of the minimum decreases 
with the system length up to the point the oscillation pattern acquires a new stripe. 
It would be interesting to experimentally test this dependence of the average MinD 
distribution on the cell length. Due to fast bleaching of the GFP we were not able to 
perform this experiment. 

In the model, the transition of MinD from one half to the other can be understood 
as follows. If MinD is localized in one half, MinE will bind and drive MinD off 
the membrane. Although the distribution of cytosolic MinD is homogenous, MinD 
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Figure 3. Oscillatory solutions of the dynamic equations (|10(1 and pill . a,b) 
Space-time plots of the total MinD and MinDE distributions on the membrane, 
Cd + Cde = {cd + Cde)/cmax and Cde = Cd e /c max , respectively, for system size L = 2/im. 
Both distributions show pole-to-pole oscillations with a temporal period of about 80s. 
c) The total MinD and the MinDE distribution averaged over one temporal period 
shown in (a) and (b), {cd + Cde) and {cde}- Both distribution display a clear minimum 
at x = Lq/2. d) Space-time plot of the total MinD distribution on the membrane, 
Cd + Cde for system size 2Lq. The pattern has doubled as compared to the pattern in 
the system of length Lq. Parameters are k\ — 2.1 • 10 _6 /xm 4 s _1 , &2 = 2.5T0~ 7 ^m 6 s -1 , 
and the remaining values as in Fig. [2h- 

preferentially binds in the other half, because there are more available binding sites. 
The resulting inhomogeneity of membrane-bound MinD is then amplified by MinD 
aggregation. As a consequence of the homogenous distribution of cytosolic MinE, the 
spatial dependence of the attachment rate of MinE follows the profile of membrane- 
bound MinD, and the distribution of MinDE complexes is similar to the one of MinD 
on the membrane, see Fig. EK, b. In particular, the positions of the maxima of Q e are 
linked to the position of the maxima of Cd- In the example given in Fig. EK, b, maxima 
are always located at the boundaries x = and x = L. 

As the system size is increased, the patterns change and striped patterns for q and 
Cde appear, see Fig. Eli. This reflects the finite wave number of the critical mode. In 
addition to changes in the oscillation pattern, the temporal period, too, changes as the 
system size is varied. It increases monotonically with the system size, but at certain 
sizes jumps back towards a lower value, see Fig. E^,. The discontinuities occur for the 
system sizes where the oscillatory pattern acquires a new "stripe" . For the parameter 
values used here, a new stripe appears for a system size of 3/im. 
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Figure 4. Oscillations of MinD-GFP in E. coli. a-d) Fluorescence images of MinD- 
GFP in a cell at subsequent time points separated by 20s. e) Time-average of all 
frames during one oscillation period. Two subsequent frames are separated by Is. f) 
Fluorescence intensity I obtained from a line scan of the fluorescence signal in (e). 
The background signal has been subtracted from the total signal which has then been 
rescaled with the maximum intensity during the oscillation. The slight asymmetry is 
due to bleaching during the observation period. Scale bar: 1/im. The cell length is 
L c = 2.3/im. 



We measured the temporal period of the oscillations in E. coli containing MinD- 
GFP, Fig. (see Methods). The periods fall in the range of 50s to 120s, even for 
bacteria of 10/im in length. The data indicate large variations of the oscillation period 
for cells of approximately the same length. This might be due to differences in the 
MinD and MinE concentrations for different bacteria and thus reflect individuality of 
the cells. An experimental verification would require the measurement of the protein 
concentration in an individual cell together with the temporal period of the Min 
oscillations. The oscillation periods found for the dynamic equations (fTUJ) and (fTTj) 
span the same range as the experimentally observed. Furthermore, experimentally we 
observed striped oscillation patterns only for bacteria longer than 3/xm, however, there 
is no sharp transition length for which the pattern changes. This behaviour, too, could 
be due to variations in the protein densities between different bacteria. 

In the model, the temporal oscillation period also depends on the total MinD and 
MinE concentrations, T> and £, see Fig. Et, d. It increases monotonically with the 
amount of MinD until it starts to descend slightly. As a function of the number of MinE 
molecules, the period decreases. Both dependencies are compatible with experimental 
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Figure 5. Dependence of the oscillation period T on parameters, a) Oscillation period 
of solutions to the Eqs. I|1U|) and (|11|) as a function of the system length. Black dots: 
oscillation pattern as in Fig. red crosses: oscillation pattern as in Fig. [3p, green 
and blue dots: oscillation pattern with three and four stripes, respectively. For the 
system length where the oscillation pattern changes, the period shows a discontinuity, 
b) Oscillation period measured for E. coli. The different symbols refer to the same 
oscillation patterns as in (a), error bars are of about the size of the symbols, c) 
Oscillation period in the model as a function of the average total MinD concentration 
T>. The period initially increases and then decreases slightly with T>. d) The same as 
(c) but for the average total MinE concentration £. The period decreases with the 
amount of MinE. For (a), (c), and (d), the parameters values are as in Fig. |Hk> the 
system length in (c) and (d) is 2fim. 



observations where the period has been found to increase with the MinD concentration 
and to decrease with the MinE concentration ^U|, but only few data points have been 
reported and further measurements are necessary. 

In conclusion, the solutions to the dynamic equations presented here are compatible 
with the experimental data, but further experiments are necessary in order to verify the 
discontinuous dependence of the oscillation period on the system length. 

4. The MinE-ring 

Early experiments indicated an accumulation of MinE close to the cell center [20J. 
This accumulation is commonly referred to as the MinE-ring. More recently, using 
deconvolution techniques, MinE has been found to be arranged in a helix with 
accumulation close to the cell center and, although weaker, at the cell poles [Hj. It 
has been suggested that the helical arrangement of MinE is induced by the helical 
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arrangement of MinD and that the accumulation of MinE occurs at the ends of the 
MinD helix In cells mutant for MinE, oscillations have been observed in the 

absence of a MinE-ring [21]. In that case, the temporal period is larger than in non- 
mutant cells. Still, this experiment clearly shows that the MinE-ring is not necessary 
for the oscillations. 

In the one-dimensional description presented above, MinE-rings correspond to 
maxima in the MinDE distribution. In the examples given so far, such maxima only 
occur at the system boundaries. For system lengths close to the value at which the 
pattern acquires a new stripe, maxima can be detected closer to the system's center. 
However, this is unlikely to be the mechanism for MinE-ring formation in E. coli, 
because no dependence of the existence of the ring on the cell size has been reported. 
Furthermore, as argued above, in the limit of homogenous cytosolic MinD- and MinE- 
distributions, maxima in the MinDE-distribution are induced by maxima in the MinD- 
distribution. MinD-rings are not observed experimentally, though. 

There are at least three other possible mechanisms that can in principle account 
for the observed accumulation of MinE at the ends of the MinD-helix. In the first 
mechanism, the diffusion length of cytosolic MinE, Ie = (E ) E/uEC maiX ) 1 / 2 , is shorter 
than half of the cell length. In this case, cytosolic MinE will predominantly attach 
before it has reached the opposite cell pole, which might lead to an accumulation close 
to the cell center. This seems to be the mechanism of MinE-ring formation in the models 
proposed in Refs. P21EB]- Accordingly, the ring vanished in JH] when the attachment 
rate of MinE was reduced, leading to an increase of Ie- To test whether this mechanism 
is supported by the Eqs. (H}-®, we studied the system for finite values of D D and 
De- In this case the cytosolic distributions cd and ce are not homogenous and all four 
equations have to be solved simultaneously. For the parameter values considered above, 
the oscillation patterns do not change significantly as long as the diffusion constants 
Dp and De are larger than 0.1/xm 2 /s and no maxima of Cd e independent of maxima 
of q were found. The diffusion length l E is also influenced by the value of oje- For 
Do = De = 2.5/im 2 /s, the values expected for diffusion in the cytosol, and values of 
uje smaller than 3.2 ■ 10 _4 /im/s the same behavior was found. Still larger values of u>e 
destroy the oscillations. Note, that by assumption, the one-dimensional description is 
only appropriate if the diffusion length Ie is larger than the cell diameter, i.e., Ie > 1/im. 
We conclude, that this mechanism is not supported by the dynamic equations presented 
above and can be tested only in a three-dimensional description. 

Two other mechanisms of minE-ring formation are suggested by studies of kinesin- 
subfamily Kinl3 members j22j |2H]. These proteins induce the depolymerization of 
microtubules. In this process they accumulate at both ends of the microtubule. As MinE 
might act on MinD filaments in much the same way, accumulation of MinE could follow 
from a similar mechanism as accumulation of the Kinl3-kinesins. The latter could be a 
consequence of a higher affinity of the microtubule end for binding the motor. Related 
ideas for the binding of MinE to MinD have been proposed in ^H] and also in |17j . 
The analogy with Kinl3-kinesins offers still another explanation for the accumluation 
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of MinE, namely dynamic accumulation due to processive depolymerization j21]. The 
present framework for studying the dynamics of Min-proteins is not suited for studying 
these effects, as filaments are not explicitly incorporated. Work in this direction is in 
progress. 



5. Conclusion and outlook 



We have presented a phenomenological description of the dynamics of MinD and 
MinE in E. coli. The description is based on the binding of MinD to the cytosolic 
membrane, recruitment of MinE to the membrane by membrane-bound MinD, MinE- 
induced detachment of MinD, as well as an interaction between molecules bound to 
the membrane. For a sufficiently strong attraction between membrane-bound MinD- 
molecules, these processes generate pole-to-pole oscillations of the Min-proteins. The 
phenomenological form of the current for membrane-bound MinD used in the present 
work captures generic features of the protein interaction and does not refer to a specific 
microscopic mechanism. It allows for a quantitative comparison between the oscillatory 
solutions of the dynamic equations and experimental findings. In agreement with the 
latter, oscillations with a temporal period from 40s to 120s can be obtained. This value 
is essentially determined by the detachment rate Ud e of MinDE-complexes. For the 
parameter values given in the text, the oscillatory pattern acquires a second stripe for a 
system size of 3/im, which agrees well with the smallest bacterial length for which period 
doubling is seen in Fig. |3Jd. This length is essentially determined by the ratio of the 
parameters k\ and &2- Finally, the time-averaged MinD-distribution shows a minimum 
at the cell center as we observed experimentally. Starting from an almost homogenous 
average distribution, the depth of the minimum increases with the system length. In 
E. coli, this feature could obviously be used to couple assembly of the Z-ring to the cell 
length and hence to control the cell cycle. 

The phenomenological description of the aggregation current can be related 
to microscopic descriptions of the protein dynamics. A simple process leading to 
aggregation is based on short-range pair interaction potentials. In this case, the 
phenomenological parameters are linked to microscopic quantities by 



C max k BT 

k 2 = 1 } Dd Ur 2 (17) 

and analogously for k\ and f?2- Here, U measures the strength of the MinD- interaction 
potential, r is a typical length scale of the interaction, and ksT thermal energy. 
These relations are valid whenever r is much smaller than the diffusion length Id = 



Dd/ujECmax. Assuming, as done above, a diffusion constant of 0.06/im 2 /s for membrane 
bound MinD, which falls well in the regime of measured diffusion constants for membrane 
proteins the values of the phenomenological coefficients used above imply values of 
for the interaction strength between membrane-bound MinD and 20ksT between 
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MinD and MinDE complexes. The range for MinD-MinD interactions is then 350nm 
and for MinD-MinDE interactions lOnm. While all other values are acceptable, the 
range for MinD-MinD interactions is too large for pure electrostatic interactions. This 
points to more involved microscopic dynamics of membrane-bound MinD than discussed 
here. 

Our analysis of the dynamic equations JIJ-tBJ) nas focused on the case of 
homogenous cytosolic distributions of MinD and MinE, and eg. Solutions in this 
limit are very similar to solutions to the full equations if the diffusion constants of 
both MinD and MinE have the realistic value of 2.5/im 2 /s. This implies that the 
approximation of constant cd and ce is appropriate. Apart from providing a reduced set 
of equations that is more convenient to study than the four equations of the full system, 
this approximation might also have an important implication regarding experiments. 
One might expect that oscillations should be observable in a purified system containing 
essentially only MinD, MinE, and phospholipid vesicles. The analysis presented here 
suggests that oscillations will show up in presence of a homogenous distribution of 
cytosolic proteins. Therefore, the closed geometry of the bacterium might not be 
essential and an open geometry could be used instead. A second implication of our 
analysis is that the number of available binding sites might need to be limited in order 
to produce oscillations. 

Other mechanisms that have been suggested for the Min-oscillations agree in the 
essential assumptions with the one studied here, namely the ability of ATP-dependent 
binding of MinD to the membrane, the recruitment of MinE to the membrane by MinD, 
and the release of MinD from the membrane driven by MinE. The proposed mechanisms 
differ, however, in essential points. Meinhardt and deBoer suggested that protein 
synthesis might be an essential element [T5] . which is not supported by experiments 
where the synthesis of proteins was interrupted and the oscillations still continued |10j . 
Howard et al. assume that MinD and MinE form complexes in the cytosol and bind 
together to the membrane |XB] . This implies in particular an exponential increase of the 
temporal period of the oscillations with the system length, with a period of 1000s for 
a system of length 7/zm. This is qualitatively different from the behavior reported for 
the mechanism studied in this work, see Fig. 0Ji. The experimental data presented in 
Fig. show oscillation periods that do not exceed 120s for bacteria of a length up to 
lOyum. However, more experiments are needed, in particular to obtain simultaneously 
values for the protein densities and the oscillation period of individual bacteria. 

The system studied by Huang et al. differs from the one studied here in the way 
MinD-aggregates are formed on the membrane [IB]. There MinD aggregation follows a 
one-step process: attachment to the membrane occurs with a higher rate at locations 
where MinD is already bound. In contrast we considered a two step-process, namely, 
cytosolic MinD binds first to the membrane and only then self-assembles into a filament. 
This difference might at first sight seem minor. However, it leads to striking differences in 
the model behaviours. First of all, assuming a one-step process for MinD aggregation, a 
three-dimensional geometry as well as a finite ATP-exchange rate is required to generate 
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striped oscillation patterns in long systems. Secondly, in the model by Huang et al. there 
are no oscillatory solutions at all for homogenous cytosolic distributions. Furthermore, 
as discussed above, in the model studied here, MinE-rings have not been found to form by 
the mechanism underlying formation of the MinE-ring in the model by Huang et al. The 
differences in the mechanism for MinD aggregation thus lead to striking consequences 
for the collective behaviour of the Min proteins. One possibility to discriminate between 
the two mechanisms is by studying the dynamics of Min proteins that are not confined 
to a cell. Further analysis of the models might lead to other possible key experiments. 

In combination, all proposed mechanisms underlying the Min oscillations suggest 
new experiments that will allow us to understand the Min-oscillations better. In order 
to make even closer contact with experiments, the formation of MinD helices must 
be included. Fluctuations due to the moderate number of Min-molecules might play 
an important role. First attempts in studying the influence of fluctuations on the 
oscillations have been undertaken [2Z1I2H1, but further work is needed and will probably 
yield results of relevance for pattern formation in the presence of noise beyond the 
Min-system. 

6. Glossary 

Mini-cell. DNA-free small cell that is produced by E. coli dividing close to a cell pole. 

Min proteins. Proteins involved in the determination of the division site. Mutations 

in these proteins lead to the formation of mini-cells. 
Linear stability analysis. In a linear stability analysis the stability of a stationary 

state against small perturbations is assessed by linearizing the dynamic equations 

with respect to the stationary state. 

Reaction diffusion system. Several reacting substances that are transported in space 
through diffusion. The reactions can induce instabilities of a stationary homogenous 
distribution leading to the formation of spatio-temporal patterns. 

Kinetic hopping model. Particles are confined to the sites of a lattice. Motion of the 
particles is described by hopping between sites. 
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Appendix A. Dynamics in three dimensions and reduction to one dimension 

In this Appendix it is shown, how the dynamics in three spatial dimensions can 
effectively be reduced to a description in one spatial dimension. The bacterium is 
conveniently approximated by a cylinder with radius Rq and length L. The volume 
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densities of cytosolic MinD and MinE at a given point are co(r, x) and c E (r, x), 
respectively. Here, r and i? denote the radial and azimuthal coordinate, respectively, 
while x is the coordinate along the long axis. Their time evolution is goverened by 

d t c D (r, ft, x) = -uj D (c max - c d (#, x) - c de (tf, x))c D (r, x)S(r - Rq) 

+ u} de c de ($,x)5(r - R ) + D D A 3d c D (r,$,x) (A.l) 

d t c E (r,-&,x) = -u E c d ($,x)c E (r,i9,x)5(r - R ) + u de c de ($, x)5(r - R ) 

+ D E A 3d c E (r,#,x) . (A.2) 

Here, c d and c de are the surface densities of membrane-bound MinD and MinDE- 
complexes, A 3d is the three-dimensional Laplace-operator and the factors of 5(r — R ) 
restrict attachment to and detachment from the cytoplasmic membrane to a region 
adjacent to the cell wall. 

2 

Since the diffusion constant of cytosolic MinD and MinE is of the order of 1^-, 
whereas the period of the oscillations is about lmin, it is reasonable to assume the 
density of cytosolic MinD and MinE to be homogenous perpendicular to the bacterial 
long axis. The volume densities of cytosolic MinD and MinE can then be replaced by 
surface densities c E and c E with 

c D (r, 0, x) = — c D (0, x) (A. 3) 

Kq 

c E (r,&,x) = —c E {'d,x) . (A.4) 

Kq 

Then, the equations governing the evolution of the protein densities read 

d t c D = - ^(c max -c d - c de )c D + uj de c de + D D A 2d c D (A. 5) 

Kq 

d t c E = - ^-c d c E + u de c de + D E A 2d c E (A. 6) 

Kq 

d t c d = ^(c max -c d - c de )c D - ^-c d c E - V • j d (A. 7) 

Kq Kq 

d t c de = ^-c d c E - u de c de , (A. 8) 

Kq 

where j is the aggregation current of MinD on the inner cell membrane and A 2d the 
two-dimensional Laplace operator on the cylinder surface. 

It has been shown that MinD forms a filamentous structure on the inner cell 
membrane [Hj. Projection on this structure yields line- densities, e.g., c d (x) = 
Jq 77 c d ($, x)R d-d. They are connected to the surface densities via 

cd(^,x) ^ -^—c D (x) (A.9) 
IixKq 

ce(^x)^-^—c e (x) (A.10) 
ZttKq 

c d (tf, x) « c d (x)5($ - (A.ll) 

c de ($,x) ^c de (x)5($-$(x)) , (A.12) 

where parametrizes the MinD-helix on the inner cell membrane. The dynamic 
equations for the line densities c~r>, c E , c d , and c de are then given by Eqs. (JU)-(|IJ). The 
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current jd appearing there is obtained by projection of the surface current on the 
x-direction. Note, that a description of the formation of MinD-helices would require 
also a specification of the perpendicular component of the current j^. 



Appendix B. Methods 

Bacteria of the E. coli K12 strain JS964 were generously donated by J. Lutkenhaus, 
University of Kansas. Bacteria taken from the freezer were grown overnight in 3ml 
Luria-Bertani (LB) medium at 37°C together with 3/il spectinomycin. Of the overnight 
culture 500/xl together with 50/il spectinomycin were given in 50ml LB medium and 
grown for two hours at 37°C. Expression of MinD-GFP was induced by 50/il IPTG 
and growing the bacteria at 31°C for at least one hour. Bacteria were immobilized for 
fluorescence imagery by using silane-coated cover slips. Fluorescent images were taken 
at room temperature with an inverted microscope (Axiovert 200M, Zeiss) using a CCD 
camera from Spot Diagnostic Instruments, Inc. driven by Metavue, Universal Imaging. 
The frame rate for measuring the time-average in Fig. |3] was 1Hz and varied between 
0.33Hz and 1Hz for the data in Fig. Et 1 - Data were analysed using Metamorph, Universal 
Imaging. 
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